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ABSTRACT 

An  asymptotic  analysis  is  developed,  describing 
ignition  of  a reactive  solid  suddenly  exposed  to  hot 

gases.  It  is  found  that  in  the  limit  of  high  activation 

i 

energy  of  the  Arrhenius  exothermic  reaction,  the  chemical 


heat  release  occurs  in  a thin  diffusive-reactive  layer 
near  the  surface.  The  solution  for  this  inner  layer  is 
then  matched  to  the  solution  for  the  outer  transient- 
diffusive  region.  An  explicit  formula  for  the  ignition 
time  is  obtained. 
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1.  INTRODUCTION 

Several  theoretical  analyses  of  the  ignition  mech- 
anism of  a reactive  solid  have  been  developed  in  the 
past,  (see  review  by  Merzhanov  and  Averson^).  Recently, 
asymptotic  techniques  have  been  successfully  applied 
to  describe  ignition  under  different  heating  mechanisms. 
Linan  & Williams  analyzed  ignition  of  a reactive  solid 

by  a constant  radiant  energy  flux  absorbed  at  the  sur- 
j 2 3 4 

7 face  or  in  depth  . Bush  & Williams  studied  the  influ- 

ence of  heat  losses  by  conduction  to  the  ambient  gas 
on  the  radiant  ignition  process.  Williams  and  Niiooka"* 
studied  ignition  by  forced  convective  heating  or  riew- 
tonian  heating,  and  Linan  and  Williams**  have  recently 
completed  the  analysis  of  ignition  by  a hot  plate  (sur- 
face temperature  constant) . 

In  the  present  study  we  analyze  the  ignition  proc- 
ess of  a reactive  solid  by  conductive  heating  from  a 
hot  ambient  gas,  so  as  to  complete  our  understanding 
of  ignition  of  a reactive  solid  under  all  the  main 
types  of  heating  mechanisms.  The  analysis  is  also  ap- 
plicable to  ignition  of  a reactive  solid  by  a hot  body 
when  the  thermal  responsivities  of  both  materials  are 
of  the  same  order. 
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In  addition,  this  study  is  a further  step  on  our 
analysis  of  ignition  of  a condensed  material  by  hot 
gases,  with  different  locations  of  the  exothermic  reac- 
tion. In  reference  (8)  we  assumed  an  heterogeneous 
reaction  occuring  at  the  surface  between  the  condensed 
fuel  and  the  ambient  oxidizer.  To  analyze  ignition  by 
a gas  phase  reaction,  it  is  necessary  to  analyze  first 
the  gasification  process  that  produces  gaseous  fuel 
to  react  with  the  ambient  oxidizer.  We  have  recently 
analyzed ^ this  gasification  process  and  at  the  present 
time  we  are  using  the  temperature  and  mass  fraction 
profiles  resulting  from  this  analysis  to  study  gas  phase 
ignition. 

These  analyses  have  been  developed  under  similar 
assumptions,  and  differ  only  on  the  location  of  the 
exothermic  reaction.  It  will  now  be  possible  to  use 
experimental  results  in  trying  to  ascertain  the  location 
of  the  reaction  producing  ignition  for  a given  material 
under  some  given  conditions. 

2.  FORMULATION 

We  assume  that  no  gasification  of  the  solid  occurs 


a 


during  the  ignition  transient.  Under  assumptions  such 


become 


Where  the  second  term  in  the  right  hand  side  of 


equation  (1)  represents  the  effect  of  an  exothermic 


are  defined  in  the  nomenclature 
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For  small  values  of  T the  effects  of  the  heat 
release  can  be  neglected  so  that  the  temperature 
profiles  take  their  inert  values: 


0Ig  - -61/T  - erf  (5/2/7) 


(5) 


The  analysis  that  follows  is  based  on  the  assumption 

that  the  nondimensional  activation  energy  g/T  is  a large 

parameter.  We  consider  T to  be  of  order  unity,  although 

the  results  are  also  valid  when  T is  small.  Since  g/T 

is  large,  small  increments  of  order  T/ g in  the  nondimensional 

temperature  0,  above  its  inert  value  0^,  produce  changes 

in  the  reaction  rate  by  a factor  e,  and  these  small 

increments  are  sufficient  to  produce  a thermal  runaway 

in  surface  temperature  at  a finite  ignition  time  x. 

ign 

We  anticipate  that  the  time  x required  to  produce  these 

increments  is  of  order  unity  and  therefore  x . is  of 

ign 

order  unity. 

From  equation  (5)  it  is  observed  that  the  inert 
temperature  in  the  condensed  phase  decreases  towards 
the  interior,  so  that  the  exothermic  reaction  only  occurs 
in  a small  layer  near  the  surface  where  g£/r/7  is  of 
order  unity.  Thus,  two  different  regions  exist  in  the 


condensed  phase:  an  inner  layer  with  characteristic 
space  dimension  z=35, in  which  the  reaction  takes  place, 
and  an  outer  region  with  characteristic  space  dimension  C, 
To  analyze  these  regions,  we  use  the  following  expansions 


Inner : 


,,  Fz  , ,r  , , ,r  . 2 . rz3  . 

e= + (g-)xi  +(»-)  x 2 + + • 

g/irT  3 3 1233t/ttT 


Outer : 


6=  -T  erf  (-M  + X,  + .. 

2/F  3 1 


Gas  phase: 


0 = erf  (-s— ;)  + — X_,  + ••• 

g 2/r  6 81 


The  variables  » X2  are  functions  of  z and  t and  the 

variables  X,  and  x are  functions  of  £ and  t. 

1 81 

Introducing  these  expansions  in  equations  (1) 
and  (2),  expanding  and  collecting  like  powers  of  0, 
the  following  equations  are  obtained 


r'-'r 


(9) 


3X 


81 


32X 


3 r 


£± 


(i2) 


From  the  boundary  and  initial  conditions  (3)  and 
(4)  the  following  conditions  are  obtained 


9*1 

3X2 

3z  s = 0 

3 z 

ss  — - -g  1— 

s 35 

s 

(13) 

Xj(t  ,0)=xgi  (t,0)  , Xi  (t  »°°)=Xgl  (t  ,°°)=0  (14) 


7- 


dx2  dX; 

4z  dz 


+ ^-exp  (Xi  -z//rr)- 

s ^2  L s 


'7T  T 

r2 


exp  (xls)  (17> 


*i-  - 


1 


3Xj 

HI 


exp  { -52/4 (t-t") } 


dx ' 


' T *"T 


(18) 


gl 


1 

r*T  3y 

s 

3C--J 

/t-t" 

exp  {-52/4(t-t') }d~‘ 


(19) 


vo 


From  equations  (16)  and  (17)  the  following  matching 
conditions  with  the  solution  for  the  outer  zone  are 
derived 


Xi (t  ,0)=x 


1 s 


(20) 


3Xi 

3x, 

_ n ^ . 1 

— r> 

8*2 

/itT 

s 3 z | 

= 1 

00 

3 z 

s 

r 

exp  (xlg) 


(21) 


From  equations  (18),  (19)  and  using  equations 
(13),  (14),  (21),  the  value  of  the  surface  temperature 

deviation  from  its  inert  value  is  found  to  satisfy  the 
integral  equations 


Combining  these  equations  to  eliminate 


a new 


the  following  integral  equation  for  x 


Equation  (24)  has  been  numerically  integrated 
by  a procedure  paralleling  the  one  used  in  reference  (8) 


The  results  of  this  integration  are  shown  in  figure  (1) 


Runaway  in  surface  temperature  is  found  to  occur  at 


a finite  ignition  time  a.  =0.393  . In 

ign 

variables  the  ignition  time  is  given  by 


c (T  -T.)  r(i+r) 


(E/RT.) 


when  the  thermal  r esponsiv it ies  are  constant.  Since 


of  T,  should  be  correctly  computed  taking  into  account 


as  was  done  in  reference  (8) 


»'■ 
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NOMENCLATURE 


A Pre-exponential  factor 
c Specific  heat  at  constant  pressure 
E Activation  Energy 

Q Heat  released  by  the  reaction  per  unit  mass  of  reactant 
R Universal  gas  constant 
T Temperature 
t Time 

x Space  coordinate  in  the  condensed  phase 

2 Stretched  space  coordinate  z=g£ 

3 Dimensionless  activation  energy  S=E (T ^ -Tq) /RT^ 

T Ratio  of  thermal  responsivities  r = / (pic)  / (p A c) 

g 

0 Dimensionless  temperature  0=(T-T.)/(T  -T.) 

3 go  3 

A Thermal  conductivity 

£ Dimensionless  space  coordinate: 

Condensed  phase  £=  -x/pAQexp (-E/RT.)/{A(T  -T  . ) } 

3 8®  J 

Gaseous  phase  £ = 'F/QAc  exp (-E/RT . ) / { p A c(T  -T . ) } 

g 3 gggoj 

p Density 

a Reduced  time  crat/{r(l+F)} 

T Dimensionless  time  tAQexp (-E/RT ^ ) / { c (TgQ-T ^ ) } 


mm 


X Reduced  temperature  defined  in  Eq(7) 


X Reduced  temperature  defined  in  Eqs(6)  and  (8) 

x 

T Mass  coordinate  ¥-/  P dx 


SUBSCRIPTS 


1 First  order  perturbation 


2 Second  order  perturbation 


)E/RT?  as  a function  of  the  nondimensio 


